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Resume 

Nous soulignons les caracteristiques du modele des quarks et ses succes dans la descrip- 
tion du spectre et des transitions de hadrons lourds ou legers, comparant brievement aux pre- 
miers principes de QCD. Pour montrer l'utilite actuelle des techniques et intuition physique du 
modele, nous presentons un modele des quarks semi-relativiste des facteurs de forme des mesons 
qui presente l'invariance d'echelle d'Isgur-Wise, et des corrections a celle-ci. Comme exemples 
d'applications, nous considerons les facteurs de forme semi-leptoniques et des disintegrations 
non leptoniques de mesons a saveur lourde. Nous posons la question de la violation de la 
factorisation, et montrons comment le modele des quarks peut donner des indications sur ce 
probleme. 



Abstract 

We underline the general features of the constituent quark model and its success in the 
description of the spectrum and transitions of hadrons made up of light or heavy quarks, 
briefly compared to first principles of QCD. To show the present usefulness of the techniques 
and intuitive insight of the model we present a semirelativistic quark model of meson form 
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factors that exhibits Isgur-Wise scaling and corrections to it. As examples of applications, we 
consider semileptonic form factors and non-leptonic decays of heavy flavor mesons. We address 
the issue of the violation of factorization, and show how the quark model can give hints on the 
problem. 



1 Success and limits of the quark model 

Following earlier remarks 1 , we will first briefly describe the main characteristics of the Con- 
stituent Quark Model (CQM), comparing it to the first principles of QCD, that we believe to 
be the right theory of strong interactions. We underline the relations and differences between 
both approaches, in particular which aspects of the CQM are understood in QCD, and also the 
successes of the CQM that still wait for rigorous explanations. 

1) The CQM assumes a fixed number of constituents in a given hadron. 

On the contrary, in QCD the number of constituents is not fixed, since a s (Q 2 ) grows with 
the distance, and in particular the number of gluons and quarks in the hadron wave function 
depends on Q 2 . 

2) In the CQM one considers constituents masses ; for the light quarks these are of the 

1 2.79 

order m ~ 0.3 GeV (from a v = — = ' ). The constituent masses, together with a non- 

2m 2M N 

relativistic approximation, results in an approximate SU(6) symmetry for the light hadrons, 
non-strange and strange, that works reasonably well, mainly in the baryon sector. 

In QCD, one starts with current masses in the Lagrangian, masses that for the lightest 
quarks u, d are of the order of a few MeV, resulting in an approximate Chiral Symmetry, that 
is dynamically broken. The constituent mass of the quark model can be identified with the 
quark mass dynamically generated by this spontaneous chiral symmetry breaking. However, 
the approximate SU(6) symmetry of the CQM, although empirically successful, does not have a 
theoretical basis. One must recall that attempts have been made to include the phenomenon of 
dynamical symmetry breaking in quark models, including dynamically generated mass, quasi- 
Goldstone mesons and degeneracy of the vacuum. These are introduced thanks to a second 
quantized chiral-symmetric version of the relativistic potential. But these attempts have not 
yet reached the stage of a complete phenomenology comparable to the CQM. 

3) In the CQM one assumes a Schrodinger equation with a flavor-independent confining 
potential with a short distance and a long distance pieces : 

V (r) = -- + \r + C (1) 

r 

(for heavy quarkonia, k = 0.5, A ~ 0.2 GeV 2 ). Empirically, one finds that a Lorentz vector 
short distance, that in QCD would correspond to one-gluon exchange (OGE), and a Lorentz 
Scalar long distance pieces are favored, as regards spin-dependent forces. Indeed, within this 
hypothesis, one has v 2 /c 2 spin-dependent corrections to the spectrum, a short distance spin-spin 
interaction interaction, 

Tr , . oei-ce2 AT . , N 167ia s cei • oe 2 c/ N ,~s 
feW = 6^ AVVW => — 6^ W (0GE) (2) 



a L • S force 
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The ratio of level differences in quarkonium would be, in the hypothesis of pure Lorentz 
vector potential (with (3 =< — >/< >) : 



* - ( L —tz Vector) 

M( Xl )-M( Xo ) 5 8 + 4/3 5 1 ; 



(5) 



while, assuming the short distance potential to be Lorentz vector and the long distance one 
Lorentz scalar : 

4 16 - 58 4 

R= £ < - (SD Vector, LD Scalar) . (6) 

5 16-2/3 5 V ' ; K ' 

Experiment gives for the cc system R(x c ) = 0.48 while for the bb one R(xb) = 0.64, 
R(x'b) = 0.58, strongly pointing to the spin-orbit force being like for a Lorentz scalar. Moreover, 
one has (M cog ( 3 Pj), M( 1 P 1 ) are independent of V T ) ■ 



M cog ( 3 Pj) = M( 1 P 1 ) 



(7) 



that is very well verified in charmonium (Table 2), where M cog is the average mass of X o, X\i 
X 2- It is important to emphasize that there is no long range spin-spin force, as it would be if 
the long distance potential were a Lorentz vector. 

These empirical hints of the quark model are well justified in QCD. Indeed, confinement 
is deduced from lattice QCD at strong coupling 2 . More realistic calculations on the lattice even 
with dynamical fermions suggest a potential of the form (1), as we can see in Fig. I 3 . Moreover, 
Wilson loop calculations with an expansion in 1/rriQ confirm the Schrodinger picture with the 
potential (1) and the form of the S • S and L • S interactions as suggested by the quark model 4 . 

4) In the CQM the vacuum is trivial. In QCD the vacuum is highly non-trivial. One 
has three main sets of phenomena : (i) Non perturbative effects connected with confinement: 
< GG > condensate, (ii) Dynamical breaking of chiral symmetry, < qq > condensate, quasi 
Nambu-Goldstone bosons n and K. (iii) Non-perturbative effects in the U(l) sector : < GG > 
condensate, non-triviality of the vacuum due to instantons, special status of the rj' pseudoscalar 
meson, possibility of CP violation in the strong interactions. Concerning (ii), as pointed out 
above, one can formulate generalizations of the quark model with spontaneous chiral symmetry 
breaking 5 . 

5) The Quark Model is a non-relativistic model, the expansion parameter being v/c. The 
mean value of this parameter is quite different in the different bound systems : 

Table I 

Quark mass, level spacing u and internal v 2 /c 2 for the different quarkonia 



Quarkonia 


m q (GeV) 


uj (GeV) 


v 2 /c 2 


bb 


5.12 


0.48 


0.13 


cc 


1.82 


0.46 


0.27 


qq 


0.3 


0.5-0.7 


0.7 


Qq 


0.3 


0.46 


0.7 



For the light quarks the internal quark motion is relativistic. In the quark model one considers 
two types of relativistic effects : (i) Relativistic corrections due to the binding, for example in 
current matrix elements like the nucleon axial coupling 

|^ = | (1-25) (8) 

where 8 = (js), has the right sign and order of magnitude (there are also 0(a s ) radiative 
corrections), (ii) Relativistic effects due to the center-of-mass motion, that come out in hadron 
form factors for example, as we will see below. 

In QCD, there are three important limiting regimes around which one can formulate very 
fruitful systematic expansions in small parameters: (i) The hard or large momentum regime, in 
which the expansion parameter is a s (Q 2 ) ~ 1/ log(<5 2 /A 2 ). (ii) Chiral symmetry limit for light 
quarks where the expansion parameters are, grosso modo, m q /A and p M /A (small momenta). 
Chiral symmetry translates into effective light hadron chiral Lagrangians. (iii) Heavy Quark 
limit where the expansion parameter is essentially A/toq. This leads to the useful Heavy Quark 
Symmetry. 

As we will see below, the relativistic effects due to the center-of-mass motion in the CQM 
make the link with the Heavy Quark Effective Theory of QCD. 

6) The CQM gives a good overall description of spectra and transitions of light and heavy 
hadrons, for the ground state and the excitations. We show a few examples in Fig. 2 (energy 
levels of first excited light baryons) , Table 2 (energy levels of heavy quarkonia cc, bb) , Table 3 
(radiative transitions of baryon isobars) and Table 4 (El transitions in the bb system). However, 
one must say that for the quasi- Goldstone bosons, mainly the pion, and the U(l) sector, the 
quark model finds difficulties. 

The quark model provides an intuitive understanding of the phenomena in terms of bound 
states and its wave functions, that is used almost unconsciously by everybody working on 
phenomena involving hadrons. 

The rigorous methods of QCD, like lattice QCD, confirm many results of the quark 
model for the ground state hadrons, and in addition can treat the pion as a Goldstone boson. 
However, these methods cannot yet give an overall view of the spectra and transitions of 
hadrons as the CQM does, especially because it handles very easily the excited states, which 
are hardly accessible to the fundamental methods. Moreover, the success of the non-relativistic 
quark model for light quarks, even if it is possibly amended by relativistic corrections, is not 
understood in terms of QCD. 

Table 2 

Quarkonia cc, bb levels (NR : non-relativistic model ; S • S : hyperfine splitting) ; * : input 
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Table 3 

Radiative transitions N* -> iV 7 (data from PDG 1994). 
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Table 4 

El radiative transitions in the bb system. 
The ratio T/k 3 oc (2 J + 1), predicted by the Non-relativistic Quark Model, is very well 

satisfied by the data. 
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2 Weak transitions of heavy mesons 

Weak decays of mesons are governed by the Cabibbo-Kobayashi-Maskawa (CKM) matrix : 



(u, c, i)7 M (l-7 5 ) 



K s K 6 
V V td Vts v tb 




that, in the Wolfenstein parametrization and expansion in powers of the Cabibbo angle A 



V us = sin^c, writes : 




V= -A 1-f AA 2 | . (10) 

V AA 3 (1 -p-irj) -A\ 2 

In this expansion, unitarity is approximate, A is of order 0(1), and the complex number 
p — ii] = e~ tS is responsible of the CP violation, with p and r\ of order 0(1). 

The different decays depend on different CKM matrix elements. For example, the semilep- 
tonic decays B — > D(D*)lv, B — > ir(p)lv, D — > K(K*)£u, K — > 7r£z/ depend respectively on 
Kb, Kfc, K s and V us . 

2.1 Heavy Quark Symmetry 

Assume a 0,g system made of a heavy quark O, and a light antiquark g. As we have seen 
above, in the Quark Model, the system is described by a Hamiltonian with a spin independent 
potential plus spin dependent terms that are inversely proportional to the quark masses : 

u P 2 , \r( n , !67ra s as! • oe 2 m Q m 

H= — + V(r) + — — 5{r) + --- p = — — . (11) 

2p 6 brnQin itlq + m 

In the infinitely heavy quark limit, for itlq — > oo, H — > finite limit : the reduced mass p — > m, 
the wave functions, binding energies, become independent of the flavor and spin of the heavy 
quark. The dynamics depends only then on the light quark degrees of freedom : H(r, p, oe 2 , m). 
This has a number of interesting consequences for hyperfine and fine splittings : 



S • S splittings : M B * = M B M D * = M D 

L • S splittings : M B ** - M B = M D ** - M D (12) 

that are qualitatively satisfied by experiment, with small corrections to the infinitely heavy 
quark limit. 

One can generalize to QCD these simple considerations, with an exact treatment of light 
degrees of freedom, gluons and light quarks 9 . For an arbitrary four- velocity of the heavy hadron 
v ti = jM ( v 2 = 1), one can write the four- momentum of the heavy quark in the form 

p% = m Q v» + W (13) 

where k is a typical residual momentum due to the binding within the heavy-light hadron wave 
function. Then, the heavy quark propagator can be approximated by 

m Q j)+k + m Q ^ jj + 1 
(m Q v + kf -m 2 Q 2vk 

The interesting feature of this expression is that the heavy propagator is independent of the 
heavy mass and depends only on the hadron four-velocity and a residual momentum k, very 
small compared to itiqv. 

One has a generalization of an atomic picture. Like in the hydrogen atom, if the hadron 
is at rest, the heavy quark acts as a static source of color (up to 1/mg corrections). However, 
here the picture is fully relativistic, since it is generalized to any hadron four-velocity, and this 
will have interesting consequences for the hadron form factors. One has an independence of the 
spin and flavor of the heavy quark, leading to a symmetry [SU(2Nf)] v , the subindex v meaning 
that there is a symmetry for any given v M '. 



This Heavy Quark Symmetry, discovered by Isgur and Wise 9 , has important consequences 
for the semileptonic form factors in the heavy-to-heavy case (like B — > D(D*)iv) and also, 
although weaker results, for the heavy-to-light case (like B — > 7i(p)£u) 10 . 

Let us define the form factors : 

/ M 2 — M 2 \ M 2 - M 2 

< V^\Pi >= (Mi + Mj) A,(q 2 ) [el - ^Lq^j - 

r* ■ a ( M 2 — M 2 \ f* ■ a 

< VMPj >= W Vjtfi e" . (15) 

The Heavy Quark Symmetry implies scaling laws for the heavy-to-light form factors 10 , that 
apply for example to the form factors B — > ix(p)iv and also to B — > K(K*), that could be 
related to the non-leptonic decays B — > K(K*)tp within the factorization approximation. The 
scaling law applies at fixed three-momentum q, small compared to the heavy quark mass, i.e. 
four momentum transfer close to its maximum value : 

<K(K*),q\J^P Q >=Cte. (\q\,M K ,M K , < m Q ) . (16) 



This implies, for the different form factors : 
/ (q),A 1 (q) 



f + (q),V(q),A (q),A 2 (q)~^ (h\ <^ M Q , q 2 = q 2 max ) . (17) 



In the heavy-to-heavy case, like in B — > D(D*)£u, stronger relations are valid for any value of 
q 2 (up to corrections in 1/rriQ, in practice l/m c ) 9 : 



(Mb + Mr>*) ^ J (M B + M D ,) (M B + M D * 



(M B +M D ,) 2 



(V4M B M D 2 ^4M B M D f (q 2 ) / 4M B M D * 2 _ 

(M B + M D ) U[q) - {Mb + M D ) [l - ~ V (M B + Afo.) ^ } ~ ^ ' " } 

(18) 

where £(i> V) is the so-called Isgur- Wise function that depends only on the product of the initial 
and final four- velocities. It satisfies the normalization £(1) = 1 (i.e. at q = or q 2 = q 2 mao ). 

These relations and the normalization are very important because they allow an almost 
model-independent determination of the CKM matrix element V c b from the decay B — > D*£u. 

The differential decay rate B — > D*£u writes 

rIV C 2 

h = i^ F (Mn*,M B ,w)\V cb \ 2 V 2 A ?(«,) (19) 



Ai | Al o 

where w is the variable w = vd* • vb = 2M b m d * — ' F(Md*, Mb, w) is a known function, t]a 
is a known short distance QCD correction, and the function £(w) differs from the exact heavy 
quark limit function £(w) by 0(1 /mo) corrections. In particular, at the normalization point, 
there are quadratic corrections 11 : £(1) = 1 + 5i/ m 2. From a linear fit £(w) = 1 — p 2 (w — 1) one 
finds (CLEO experiment 12 ) 

\V cb \ = 0.038 ±0.006 ±0.004 . (20) 



2.2 Quark model of semileptonic form factors 



This is a weak binding model that takes into account the relativistic center-of-mass motion of 
the hadrons 13 , as refered to in 5) of section 1. The model provides an intuitive link between 
the Quark Model and the Heavy Quark Symmetry. 
The total hadron wave function writes 



*p({ft}) 



*(5>-p)*p({Pi» 



(21) 



with the internal wave function 



*p({Pi}) — N 



ns(p) 



^p=o({p,}) 



The model accounts for two important relativistic effects : 
1) Lorentz contractions of the wave function 

PiT = PiT 

E 



E 



P 



M M 
2) Lorentz boost of Dirac spinors 



Viz 



ffi-P ' 
E + M 



(22) 



(23) 



(24) 



The calculation of current matrix elements is relatively simple in the equal velocity frame 
(EVF), where the initial and final hadrons have the same modulus of the velocity (3, and 
opposite in direction. An interesting relation in this frame is : 
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The form factors read, in this model : 
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(26) 



etc. In these expressions I(q 2 ) is the wave function overlap, and X + , Xy, Xi, X 2 (Mi, Mf, q 2 ) 
~ 0(m) are corrections proportional to the spectator quark mass. 

The model has a number of interesting features. It satisfies heavy-to-light Isgur-Wise 
scaling (16), (17), and also heavy-to-heavy Isgur-Wise scaling (18) with the Isgur-Wise function 



av ■ v') 



-I(vv') £(1) = 1 



1 + v ■ v' 

where I(v ■ v') depends on the overlap of the wave functions at rest : 

m \ _ / m 



I(v ■ v') = J $J /=0 
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where q is the Lorentz contracted momentum transfer 
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(27) 



(28) 



(l-^W = (Mf + MiY 



f-q 2 
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2 / V ■ V 



V • V' + 1 



(29) 



The normalization £(1) = 1 results then from the normalization of the wave function at rest. 
In the harmonic oscillator model one has : 



av ■ v f ) 



1 + V • V 



7 exp 



m 2 R 2 (v-v' — l 



V2 \v-v' + l 



an expression found elsewhere on different grounds 14 . The slope is given by 

2 1 m 2 R 2 



(30) 



(31) 



where the second term is the dominant one in the non-relativistic limit, and the first is a 
relativistic correction due to the boost of spinors. 

Moreover, the model exhibits scaling corrections due to the spectator quark mass that 
point to a softening of the scaling laws (17) at fixed q. These softening of the scaling is confirmed 
by lattice calculations within very large errors 15 . Also, neglecting hyperfine splitting, the model 
gives relations among the form factors that express the fact that, in the model, the total quark 
spin of the final meson equals the meson spin : 



A (q 2 ) = f + (q 2 ) 



2M f (Mi - M f ) f (q 2 ) = (M 2 - M 2 - q 2 ) A x (q 2 ) 



(32) 



\(M 2 ,M 2 ,q' 
(Mi + M f ) 2 



A 2 (q 2 ) ■ 



2.3 Example of phenomenological application 



Let us consider the relation between the semileptonic form factors D — > K(K*)£u and the 
form factors extracted from the non-leptonic decays B — > ipK{K*). We have the following 
theoretical constraints and experimental data : 

1) Heavy-to-light Isgur-Wise scaling near q 2 = q^ax- 

2) Data onD^ K(K*)£u near q 2 = 16 : 



/; c (0) = 0.77 ±0.04 
Af (0) = 0.61 ±0.05 



V sc {0) = 1.16 ±0.16 
A S 2 C (0) = 0.45 ± 0.09 



3) Data on the decays B — > ^/K(K*), that give information on the form factors at q 2 



(L stands for longitudinal polarization) : 

R _ r (gj -» ^*°) 
r (s° - 



Ri 



(33) 



(34) 



Table 5 





R 


Rl 






0.66 ±0.10 + ° n l° Q (CDF) 21 
— U.Uo 


Experiment 


1.64 ± 0.34 (CLEO) 23 


> 0.78 (ARGUS) 22 
0.80 ± 0.08 ± 0.05 (CLEO) 22 


BSWI 17 


4.23 


0.57 


BSWII 18 


1.61 


0.36 


GISW 19 


1.71 


0.06 


QCDSR 20 


7.60 


0.36 



Within the factorization hypothesis 17 these rates can be related to the B 
factors, like for example : 

G 



A - j,K) = ^V cb V; s 2U m B f + (ml)a 2 p 



K(K*) form 
(35) 



Table 6 





R 


Rl 


/• c (0) 
A\° (0) 


V sc (0) 
Af c (0) 


Af(0) 
A\ c (0) 


X 2 /DoF 


Exp. 


1.64 ± 0.34 


0.66 ± 0.14 


1.26 ± 0.12 


1.90 ± 0.25 


0.74 ± 0.15 




Fit 


2.15 


0.45 


1.45 


1.62 


0.81 


4.2 



A simultaneous fit to D semileptonic form factors and to B non-leptonic decays is not very 
good (Table 6), but some trends of our model seems to be favored, namely : 



fsb( 2 \ 
J \HmaxJ 

J 4- \1max J 



rn B 



-l»K 



m D +m K ) ' 



etc. 



1) Softened scaling : 

2) Pole for the ratios f+(q 2 )/A 1 (q 2 ), V \q 2 ) / 'A^q 2 ) , A 2 (q 2 ) / A 1 (q 2 ) 

3) Weak q 2 dependence of A 1 (q 2 ). 



2.4 The factorization problem 



In the precedent application, we have made use of factorization to establish a connection be- 
tween non-leptonic decays and semi-leptonic form factors. For B mesons, factorization seems 
approximately correct at least for the so-called class I decays, that are of leading order at large 
N c , as shown by CLEO results 3 . But it could happen that for the decays B — > tpK(K*), that 
are subleading in 1/N C , violations of factorization (i.e. the recipe (35)) could modify the ratio 
R and even R^. We can wonder whether this assumption has theoretical justification. First, 
we recall that this factorization principle is implemented in the standard approach of SVZ 24 , 
with particular values of a\ and a 2 , which are found to be 



ax — ci + 



C-2 



a 2 = c 2 + 



ci 



where c±, c 2 are the coefficients in the effective weak Hamiltonian : 



(36) 



H = c x {cb)i{sc)i + c 2 (cc) 1 (sb) ] 



(37) 



In fact it can be shown that this standard factorization is exact in the N c —>■ 00 limit. However, 
as shown by a simple argument due to Shifman and by us 25 , the predicted rates cannot be 
correct at subleading order 1/N C , because they contradict duality. The argument is as follows. 
Let us consider as an example the Cabibbo suppressed process b — > cus. There are two types of 
decays corresponding to different color topologies, namely B d — > D + D~ (class I) and B d — > ipK° 
(class II). Let us consider the effective Hamiltonian (the subindex means color singlet). The 
decay rates into a pair of hadrons corresponding to classes I and II result, from factorization : 



\A{B S 
\A(B S 



,sc Acs h 



[cchiss 1 ~ 



l C2 

01 + K. 



C 2 + 



(31 



Total quark decay rate differs from the sum of the expressions (38), since it is proportional to : 



\A(b-,cus)\ 2 ~cl + cl + 2 C -^=£ 



Cl + 



C2_ 



+ 



C2 + 



Nr. 



1 2 



(39) 



On the contrary, duality tells us that by summing over all mesons one should obtain total 
decay rate oc c\ + cl + 2cic 2 /N c . Therefore, Shifman proposes a recipe: simply to impose 
duality by multiplying the subleading terms in the decay rates to hadrons by a factor x — 1/2 



\A(B d 



Isc 1 (cs 1 ~ 



\A(B d -> {cc)x(ss) 



Cl+X w c . 

c 2 + x- 



(40) 



In our opinion, the reason for the discrepancy observed in ( p9|) is that the color configu- 
rations (cb)i(sc)i and (cc)i(sb)i are not orthogonal, the overlap being of order 1/N C : 



(cfe)i(sc)i = — (cc)i(s6)j + 



N? - 1 



Nr 



[(cc) 8 (s6) 8 ] 1 



(41) 



We have considered a model with non-relativistic scalar quarks bound to color singlets by 
a color harmonic oscillator potential. This model takes into account of a Final State Interaction 



that automatically restores duality or the conservation of probability. Indeed, summing over 
all mesons in the limit in which the radius R — > oo one should find the free quark result. This 
is indeed the case, since one finds : 



{AiB^iccUss)^ 2 - c 2 + (l-y) 



(42) 



with y(m c ,m s ) depending non-trivially on the masses, unlike the universal factor 1/2 proposed 
by Shifman. It is seen that FSI restores automatically duality. Notice by the way that expression 
(42), for an arbitrary value of y, is more general than Shifman's Ansatz and would as well restore 
duality. Another study of restoration of duality in the Quark Model approach may be found in 
[26]. The ratios R and R L could be modified by this type of FSI effects, since by introducing 
spin there is no reason a priori that the ratio Rl would not be modified also. 
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